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Abstract 

We perform a detailed study of the renormalization group equations in the inverse seesaw model. 
Especially, we derive compact analytical formulas for the running of the neutrino parameters in 
the standard model and the minimal super symmetric standard model, and illustrate that, due to 
large Yukawa coupling corrections, significant running effects on the leptonic mixing angles can 
be naturally obtained in the proximity of the electroweak scale, perhaps even within the reach of 
the LHC. In general, if the mass spectrum of the light neutrinos is nearly degenerate, the running 
effects are enhanced to experimentally accessible levels, well suitable for the investigation of the 
underlying dynamics behind the neutrino mass generation and the lepton flavor structure. In 
addition, the effects of the seesaw thresholds are discussed, and a brief comparison to other seesaw 
models is carried out. 



*Electronic address: johbergs@kth.sc 
^Electronic address: malinsky@kth.se 
^Electronic address: tommy@theophys.kth.se 



3 Electronic address: zhanghc@kth.se 



I. INTRODUCTION 



Experimental progress on neutrino masses and leptonic mixing has opened up a new 
window in searching for new physics beyond the standard model (SM) of particle physics 
during the past decade. Since neutrinos are massless particles in the SM, one usually extends 
the SM particle content in order to accommodate massive neutrinos. Among various theories 

n 

of this kind, the seesaw mechanism [1[ attracts a lot of attention in virtue of its naturalness 
and simplicity. For instance, in the conventional type-I seesaw model, three right-handed 
neutrinos with a Majorana mass term far above the electroweak scale are introduced. The 
masses of the light neutrinos are then strongly suppressed with respect to the masses of the 
charged leptons by the ratio between the electroweak scale and the mass scale of the heavy 
right-handed neutrinos. Usually, the neutrino parameters are measured in low-energy scale 
experiments, while on the other hand, the seesaw-induced neutrino mass operator often 
emerges at some very high-energy scale. Therefore, neutrino masses and leptonic mixing 
parameters are subject to radiative corrections, i.e., they are modified by the renormalization 
group (RG) running. In principle, the RG running effects can even be physically relevant, 
especially if the seesaw scale turns out to be not extremely high. 

The generic features of the RG running of neutrino parameters have been investigated 
intensively in the literature. Typically, at energy scales lower than the seesaw threshold, i.e., 
the mass scale of the heavy seesaw particles, the RG running behavior of neutrino masses 
and leptonic mixing can be described within an effective theory, which is essentially the same 
for various seesaw models, reflecting the uniqueness of the dimension-five Majorana mass 
operator in the SM. However, at energy scales higher than the seesaw threshold, a full theory 
has to be considered, and the interplay between the light and heavy sectors could make the 
RG running effects particularly different compared to those in the effective theory. The RG 
running effects above the seesaw scale are especially relevant in the wide class of popular 
theories based on the idea of grand unification, i.e. grand unification theory (GUT), where 
the specific flavor structure stipulated at the GUT scale, typically of the order of 10 16 GeV, 
often experiences a long range running over many intermediate energy-scale thresholds. 

Concerning any specific seesaw model, a generic analysis of the RG running of neutrino 
parameters is clearly inconceivable because of an infinite number of possible underlying 
models, in particular, above the relevant seesaw threshold. Nevertheless, the universal fea- 
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tures, like the presence of extra degrees of freedom underpinning the seesaw, i.e., the heavy 
fermions in the type-I and -III models or the scalar triplets in the type-II seesaw model, 
and the nature of their interactions with the SM sector, still admit for a high degree of 
theoretical scrutiny, providing a basis for any further work within a specific class of unified 
models. 

The full set of renormalization group equations (RGEs) in the type-I, -II, and -III seesaw 
models have been derived, both in the SM and in the minimal supersymmetric standard 
model (MSSM) 2|-[5|. However, in these models, the seesaw scale is often taken not too 
far below the GUT scale, which hinders the direct experimental testability of the origin of 
neutrino masses. Although, in the type-I seesaw model, one can bring down the mass scale of 
the right-handed neutrinos by means of a severe fine-tuning among the model parameters , 
radiative corrections induced by right-handed neutrinos tend to spoil the stability of such 
settings beyond the tree-level approximation. 

Recently, a lot of attention has been focused on low-scale seesaw models, and especially, 
the possibility of searching for seesaw particles at the Large Hadron Collider (LHC), see, 
e.g., Ref. {^J and references therein. In these models, the smallness of the neutrino masses 
is protected by other means than the GUT-scale suppression, such as a small amount of 
lepton number breaking. A very attractive example along this direction is the inverse seesaw 
model js], in which the light neutrino masses are driven by a tiny Majorana mass insertion 
in the heavy Dirac neutrino mass matrix instead of the mass scale of the heavy neutrinos. 
Furthermore, due to the pseudo-Dirac feature of the heavy singlets, the model does not suffer 
from either large radiative corrections or unnatural fine-tuning problems. This admits for 
bringing the heavy neutrinos down to the LHC energy range while retaining potentially large 
Yukawa couplings. Subsequently, this makes the model phenomenologically very attractive 
from the lepton flavor violation point of view (9] or for the potentially large nonunitarity 
effects in the leptonic mixing matrix [10 |. Since the RG running, and in particular, the 
threshold effects, in such a scenario can play an important role even at a relatively low 
energy scale, there is a need to look in detail at the running effects on neutrino parameters. 

In this work, we will investigate in detail the RG evolution of neutrino masses and leptonic 
mixing parameters in the inverse seesaw model. In particular, in Sec. [Til we first review 
briefly the inverse seesaw model. Next, in Sec. IIIH we present the full set of RGEs for the 
neutrino parameters. Then, in Sec. HV] we provide useful approximations of the RGEs found 
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in Sec. IIHI Section |V] is devoted to detailed numerical illustrations and interpretation of the 
RG running behavior of the light neutrino mass and leptonic mixing parameters. In Sec. IVIl a 
discussion and comparisons between RGEs in various seesaw models are performed. Finally, 
in Sec. IVII[ a summary is given and our conclusions are presented. In addition, in App. [AJ 
the complete one-loop RGEs for some of the neutrino parameters are listed. 



II. THE INVERSE SEESAW MODEL 



The inverse seesaw model is constructed by extending the SM particle content with three 
right-handed neutrinos v R = (i/ri, z/r,2, and three SM gauge singlets S = (Si, S2, S3). 1 
The Lagrangian of the lepton masses is then arranged so that it reads in flavor basis 

-C v = h<pY e e R + hA>Y v v R + S~ c M R is R + ±&M S S + h.c. , (1) 

where <fi is the SM Higgs fields with <p = \T2<\>* , M R is an arbitrary 3x3 matrix, and M5 is 
a 3 x 3 complex symmetric matrix. Here Y e and Y v are the corresponding Yukawa coupling 
matrices, and in general, they are arbitrary. Without loss of generality, due to the freedom 
in the right-handed fields, one can always perform basis transformations to keep both M R 
and Y e diagonal, i.e., M R = diag(M 1; M 2 , M 3 ) and Y e = diag(y e , y^, y T ). Since the new extra 
fermions are SM gauge singlets, their masses are not protected by the Higgs mechanism, 
and could be related with some new physics beyond the SM. Thus, one expects the scale of 
heavy seesaw particles to be higher than the electroweak scale, i.e., M R > (9(100) GeV. At 
energy scales lower than M R , heavy degrees of freedom in the theory should be integrated 
out, which leaves a series of higher-dimensional nonrenormalizable operators in the effective 
theory. At dimension-five level, the only allowed operator is the so-called Weinberg operator, 
coupling two lepton doublets to the SM Higgs field, on the form 

£f 5 = -^(40)-M0 T ^)+h.c., (2) 



1 Note that in order to accommodate the experimentally measured neutrino mass-squared differences and 
mixing angles, the minimal setup is to introduce only two right-handed neutrinos together with two heavy 



singlets [llj, which leads to one massless neutrino in the model. In this work, we will concentrate on 
the more general case with three right-handed neutrinos together with three heavy singlets, which is the 
minimal case allowing three massive light neutrinos. 
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where k is the effective coupling matrix. In the inverse seesaw model and at tree-level, it is 
given by 

K~Y v M^M s {Ml)- x Yj. (3) 

After electroweak symmetry breaking, the dimension-five operator defined in Eq. (j2J) yields 
an effective Majorana mass term for the light neutrinos 

m u = kv 2 , (4) 

with v ~ 174 GeV being the vacuum expectation value of the Higgs doublet. 

Comparing Eq. flS]) with the typical type-I seesaw model, in which the effective coupling 
matrix k is given by k ~ Y u M^ 1 Yj , one can observe that, in the inverse seesaw, the masses of 
the light neutrinos are not only suppressed by Mr, but also by the small Majorana insertion 
Ms- In the limit Ms — > 0, k = indicates the restoration of lepton number conservation, 
which in turn ensures the naturalness of a tiny Ms- 

Now, since the energy scale of Mr is not necessarily high, we may naturally ask if it is 
related to some new physics around the TeV scale and could be tested at the current hadron 
colliders (e.g., LHC). In fact, for the purpose of searching for heavy singlets at colliders, 
it is very important to have visible mixing effects between the light and heavy neutrinos. 
Furthermore, the admixture results in nonunitarity effects in neutrino flavor transitions, 



especially for some of the future long-baseline neutrino oscillation experiments |1CH14|. 

In spite of the underlying physics responsible for large neutrino masses, the particle 
content of the inverse seesaw model is essentially the same as that of the type-I seesaw model, 
but with six right-handed neutrinos. In principle, one may treat the heavy singlets Si (i = 
1,2,3) as three extra right-handed neutrinos, possessing vanishing Yukawa couplings with 
the lepton doublets. In the limit Ms — > 0, lepton number conservation is restored, and the 
six heavy singlets can be combined together to form three four-component Dirac particles, 
while keeping the light neutrinos massless. Therefore, any phenomena in association with 
lepton number violation should be proportional to the Majorana mass insertion Ms- 



III. RGES FOR THE NEUTRINO MASS MATRIX 



As discussed above, the analogy between the type-I seesaw model and the inverse seesaw 
model allows us to readily write down the RGEs for the Yukawa couplings and the neutrino 
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mass matrix from the existing RGEs of the type-I seesaw model Q . Especially, the relevant 
beta functions, obtained using the minimal subtraction renormalization scheme, are given 
by 

dY 

16tt 2 /^ = (« e + C c ;X e + C v e X v ) Y e , (5) 

16^—^ = + ^ + C v Xv) y v (6) 

d/i 

IGttV^ = C R M n YjY v , (7) 

16ttV^ = , (8) 
dfx 

where /i is the renormalization scale, Y) (for / = u,d,e,u) denote the Yukawa coupling matri- 
ces with X f = YfY}. The coefficients (C e e , C* C%, C R ) stand for (3/2, -3/2, -3/2, 3/2, 1) 
in the SM and (3,1,1,3,2) in the MSSM, respectively. 2 The coefficient a u is flavor blind 
and reads 

a v = tr (3X U + 3X d + X e + X u ) - - -g\ (9) 

in the SM, and 

a v = tr (3X U + X u ) - -g\ - 3g\ (10) 

in the MSSM with gi being the gauge couplings. Then, if we make use of m v at energy scales 
above the seesaw threshold, we can derive from Eqs. ([3]) and ©-((8]) 

dm,. 



m v = 2a v m v + {C e v X e + C m X v ) m v (11) 
+m v [C e v X e + C m X u f , (12) 



with C m = 1/2 in the SM and C m = 1 in the MSSM. Here, for simplicity, we have defined 
t = l/(167r 2 )ln(/i/M z ). 

The beta functions in Eqs. (jSJ)-© are obtained in the minimal subtraction renormaliza- 
tion scheme, and thus, we use an effective theory below the seesaw thresholds, in which the 



2 Recall that this difference is, namely, due to the fact that the proper vertex corrections are absent above the 
soft supcrsymmetry-breaking scale in the MSSM. The differences in signs turn out to be very important 
for the RG behavior of the leptonic mixing angles above the seesaw threshold, see, e.g., Ref. [15| and 
references therein. 
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heavy fields are absent and their physical effects are covered by a series of higher-dimensional 
operators. To ensure that the full and the effective theories give identical predictions for 
physical quantities at low-energy scales, the parameters of the full and the effective theories 
have to be related to each other. 

In the case of the neutrino mass matrix, this means relations between the effective cou- 
pling matrix k and the parameters Y u , Mr, and Ms of the full theory. This is technically 
called matching between the full and the effective theories. In the effective theory, the Wein- 
berg operator (J2J), responsible for the masses of the light neutrinos, does not depend on the 
specific seesaw realization and its evolution equation reads 

k = a K K + (CZX e ) k + k (C e v X e f , (13) 

where 

a K = 2 tr (3X„ + 3X d + X e ) + A - 3g% , (in the SM) (14) 

a K = tr (6X U ) - ^gl - 6g 2 2 , (in the MSSM) (15) 

with A denoting the SM Higgs self-coupling constant. 

For the simplest case, if the mass spectrum of the heavy singlets is degenerate, namely 
Mi = M2 = M3 = Mo, one can simply make use of the tree- level matching condition 

k|„ = YvM^MsiMlY'YJl (16) 

at the energy scale \x = Mq. In the most general case with nondegenerate heavy singlets, 
i.e., Mi < M 2 < M 3 , the situation becomes more complicated and the heavy singlets have to 
be sequentially decoupled from the theory For instance, at energy scales between the 
nth and (n — l)th thresholds, the heavy singlets are partially integrated out, leaving only 
a 3 x (n — 1) submatrix in Y v , which is nonvanishing in the basis, where the heavy singlet 
mass matrix is diagonal. The decoupling of the nth heavy singlet leads to the appearance 
of an effective dimension- five operator similar to that in Eq. (j2J), and the effective neutrino 
mass matrix below M n is described by two parts 



m, ( , n) = v 2 



(17) 



where (n) labels the quantities relevant for the effective theory between the nth and (n — l)th 
threshold. In the SM, the RGEs for the two terms above have different coefficients for 



the gauge coupling and Higgs self-coupling contributions, which can be traced back to the 
decoupling of the right-handed neutrinos from the theory. However, such a mismatch is 
absent in the MSSM due to the supersymmetric structure of the MSSM Higgs and gauge 
sectors. 3 Therefore, this feature may result in significant RG running effects only in the SM, 
in particular, when the mass spectrum of the heavy neutrinos is quite hierarchical. 



IV. ANALYTICAL RGES FOR NEUTRINO PARAMETERS 



In order to analytically investigate the RG evolution of leptonic mixing angles, CP- 
violating phases, and neutrino masses, one can translate the full RGEs for the neutrino 
mass matrix into a system of differential equations for these parameters. The correspond- 
ing formulas have been discussed below the seesaw scale [3], as well as above the seesaw 



thresholds in the type-I [15j, -II |4|], and -III |5[ seesaw models. 

In general, there are two different strategies for solving the resulting system. In the 
top-down approach, the initial condition is specified at a certain high-energy scale, often 
motivated by the flavor structure of a specific GUT-scale scenario. This is an advantage, 
since all the necessary ingredients are fixed at the high-energy scale and the running down 
to the electroweak scale is a mere tedium. On the other hand, only some regions in the 
parameter space of the full theory would presumably lead to good fits of the low-energy 
data, which, together with an ah initio model dependency, makes this method potentially 
quite inefficient. 

In this study, we instead adopt the bottom-up approach, in which the initial condition 
for the observables of our interest is fixed at the low-energy scale, thus exploiting all the 
available experimental information from the beginning. It is also clear from the mere pa- 
rameter counting that at the matching scale one can determine the underlying theory pa- 
rameters only up to an equivalence class defined by the matching condition, and thus, the 
model-dependency problem is somewhat delayed. Nevertheless, in order to set off from the 
threshold, a representative should be chosen, which amounts to adding extra assumptions 
on the flavor structure of the underlyin g th eory at the matching point. For a more detailed 



discussion on these issues, see e.g. Ref. 



18| and references therein. 



We assume that the relevant threshold is above the soft supcrsymmetry-breaking scale. 
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In particular, one can make use of some of the qualitative features of the RGEs of the 
full theory. For instance, in the type-I seesaw model, the structure of Eqs. (jSJ)-© an d 
(ITT]) justifies the typical diagonality assumption made on Y v . Indeed, if both Y e and Y v are 
diagonal at a common energy scale, then they will retain the one-loop diagonality at any 
energy scale, since there is no nondiagonal element in the RGEs for Y e and Y v . Thus, for the 
sake of simplicity, we assume Y v = diag(y 1/ i, y V 2, y v z) in the basis in which Y e is also diagonal, 4 
which is reflected by the fact that from now on we work with just three combinations of the 
leptonic Yukawa couplings, namely 



(2/1, vh vl) = c t(y 2 e, yl, vl) + c m (y 2 ul , y 2 u2 , viz) 



(18) 



We will briefly comment on the general case with a nondiagonal Yukawa coupling matrix 
Y v later in Sec. IV Dl The leptonic mixing matrix corresponds then to the unitary matrix U 
diagonalizing 



k = C/diag(/ci, fa, kzjU 1 



(19) 



where hi (i = 1,2,3) are the eigenvalues of k. It is usually given using the standard 
parametrization 



U 



C12C13 



S12C13 



si3e 



-A 



— S 12C23 — c 12S23 s 13 el5 C12C23 — Si2S23Si3e l5 S23 c 13 



V 



S12S23 - ci2C 23 si 3 e 



16 



-C12S23 - Sl2C 2 3Sl3e 1<5 C23C13 



7 



\ 



(20) 



V 



with Cij = cos9ij and = sin^j (ij = 12, 13, 23). Here = diagi^e 1 ^ 1 , e 1 ^ 2 , e 1 ^ 3 ) denotes 
three unphysical phases, which are required for the diagonalization of an arbitrary complex 
symmetric matrix. 

For the RG running of leptonic mixing angles and CP-violating phases, in order to simplify 
the results, we define the quantities 



c, 



m, — rrij 
mi + m.j 



(21) 



4 Note that the simplifying assumption of a simultaneous diagonality of Y e , Y Vl and Mr is specific for the 
inverse seesaw model and cannot be imposed in e.g. the type-I seesaw model. 
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Inserting Eqs. ( !T9|) and (12T|) into Eq. ( TTTj) and using Eq. ( ITB"]) . we arrive after some tedious 
calculations at the RGEs for the leptonic mixing angles and the CP-violating phases in the 
current scheme. The explicit, but rather cumbersome results, are shown in App. 

In the case that the mass spectrum of the light neutrinos is nearly degenerate, i.e., 
rrii — m 2 — m 3, one can expect large enhancement factors C^ 1 , which strongly boost the 
RG running effects. In particular, for the leptonic mixing angle #12, the main RG running 
effect comes from the correction proportional to , and we obtain approximately 

#12 - j~c p . a {sncncp-a [c\ 3 yl + (s 2 23 s 2 13 - c 2 23 ) y\ - (s 2 23 - c 2 23 s\ 3 ) yf\ 
(,21 

+S23C23S13 (cs+p-a - 2s 2 12 csc p - a ) (y 2 3 - yl) } . (22) 

If further neglecting the small leptonic mixing angle #13, we arrive at 

O12 ^ ^s 12 c 12 c 2 (yf - c 2 23 y\ - s 2 23 yl) . (23) 

In addition, for the leptonic mixing angles 9 23 and 6*13, we only keep the terms, which are 
not suppressed by s 13 , and obtain approximately 

#23 - -^ 23 C23 (4j<* + c\ 2 cl) (y 2 2 - yf) , (24) 

S31 

#13 - -^si 2 ci2S2 3 C23 (c a C5+a - c p c s+p ) (yj - yf) , (25) 

<-,31 

where the approximate relation ( 31 ~ ( 32 has been used. 

Now, we turn to the analytical RGEs for the CP-violating phases. In the limit #13 — > 0, it 
is worthwhile to mention that the Dirac CP-violating phase S loses its meaning. However, it 
has been pointed out that, with the RG running, both nontrivial 6* 13 and 5 can be generated 
radiatively [ljj]. Therefore, we keep terms proportional to either the inverse power of s± 3 or 
£21 , and obtain 

8 - — ^si2C 12 s 2 3C 2 3 {cp-ss+p - c a s s+a ) (y 2 - yf) 

s 13<-,31 

+T^—s 2p -2a (yl - c 2 23 y 2 2 - s 2 23 yf) , (26) 
A21 

P - 7^—c 2 12 s 2 p- 2 a (c 2 23 y 2 + s 2 23 y 2 3 - yf) , (27) 
sl 2 s 2 p- 2 a (c 2 23 yj + s 2 23 y 2 3 - yf) . (28) 



a 



2C: 



21 
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Finally, we express the analytical RGEs for the masses of the light neutrinos as 

m x ~ 2m 1 [c\ 2 yl + (c\ z yl + s 2 2Z yl) s\ 2 ] + 2m 1 a„ , (29) 
m 2 ~ 2m 2 [s\ 2 yl + (cl 3 yj + s^yf) c\ 2 ] + 2m2a„ , (30) 
m 3 ~ 2m 3 (s^^ + cl-iVl) + 277230,, . (31) 

One can observe that the yis play the key role in the RG running of #12 . In fact, in the 
low-energy scale type-I seesaw model, y V i <C y T has to be satisfied in order to effectively 
suppress the masses of the light neutrinos, and therefore, no visible RG running effects can 
be achieved in the SM. However, in the inverse seesaw model, the y^s can be chosen to be of 
order unity, since the masses of the light neutrinos are diminished by the Majorana insertion 
Ms instead of Y u , and therefore, sizable RG running effects can be naturally expected in the 
inverse seesaw model. 

Concerning the bottom-up approach adopted in this study, let us add one more technical 
remark at this point. Upon crossing the seesaw threshold, the matching between the full 
and the effective theories can be very easily performed in a basis, where the mass matrix 
of the heavy singlets is diagonal. However, since Mr 3> Ms is well satisfied in the current 
framework, we can effectively work out the matching in a basis, where Mr is diagonal. The 
inaccuracy induced by this approximation is limited to 0(Ms/Mr), which can be safely 
ignored in our numerical calculations. 



V. NUMERICAL ANALYSIS, ILLUSTRATIONS, AND INTERPRETATIONS 

We proceed to the numerical evolution of the RGEs in order to show the representative 
RG running behavior of the neutrino parameters in the inverse seesaw model. In our compu- 
tations, we solve the one-loop RGEs exactly instead of using the approximate formulas. For 
this purpose, we adopt the values of the neutrino mass-squared differences Am?- = m 2 — m 2 



and the leptonic mixing angles from a global fit of current experimental data in Ref. [20 1 
that we assume to be given at the energy scale \x = Mz- Note that, since there is no 
compelling evidence of a nonzero #13 so far, we set #13 = throughout the numerical illus- 
trations, unless otherwise stated. We also use the values of quark and charged-lepton masses 



as well as the gauge couplings given in Ref. 21] . In the SM case, we choose the Higgs mass 



to be m H = 140 GeV, and we assume the shape of the Higgs spectrum to be driven by 
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tan /3 = 10 and tan /3 = 50 as well as = Mz in the MSSM, where the latter assumption 
is motivated by simplicity. Currently, there is no direct experimental information on the 
absolute values of the masses of the light neutrinos. However, the recent measurement on 
the cosmic microwave background finds that the sum of the masses of the light neutrinos is 
less than 0.58 eV at 95 % CL. j^j. One can estimate that this constraint can be satisfied 
if rrii < 0.15 eV (i = 1, 2, 3). In addition, the sign of Am| 2 ~ Am^ remains undetermined. 
Therefore, we have four representative choices for the mass spectrum of the light neutri- 
nos: normal and hierarchical neutrino mass spectrum (NH) with mi <C m 2 <C m 3 ; normal 
and nearly degenerate neutrino mass spectrum (ND) with mi < m 2 < m 3 ; inverted and 
hierarchical neutrino mass spectrum (IH) with m 3 <C mi <C m 2 ; and inverted and nearly 
degenerate neutrino mass spectrum (ID) with m 3 < mi < m 2 . As we have shown in both 
the NH and IH cases, there is no visible enhancement factor, since ( 2 ~\ = C31 = C32 — 1 i n 
the NH case and C21 = — C3I 1 = — C32 1 — 1 in the IH case. Thus, we will mainly work in 
the ND and ID cases in order to gain sizable enhancement factors. For example, in the ND 
case with mi = 0.1 eV, C2I 1 ~ 500 and (^j 1 ~ C-^ 1 ~ 20 hold to a good precision. As for 
the masses of the heavy singlets, it has been pointed out that, if the masses are around the 
TeV scale, one may search for their signatures at the LHC via the trilepton channels, i.e., 



pp — > £££u, in which the SM background is relatively small [23]. We thereby choose these 
masses located around the TeV scale in our numerical studies. In practice, we first consider 
the simplest situation with all masses of the heavy singlets being identical, and then come 
to the most general case with a nondegenerate mass spectrum of the heavy neutrinos. 

However, for such a low seesaw scale, one has to consider further phenomenological con- 
straints on the parameter space of the model. In particular, significant nonunitarity effects 
in the leptonic mixing can emerge in the TeV scale inverse seesaw if there are 0(1) Yukawa 
couplings in the lepton sector ll|, [l3|. This is traced back to an effective dimension-six 
operator 

£f 6 = c (40) # (^4) , (32) 

where c = (Y U M Y{ 1 ^ (Y^MpT 1 )^ at leading order. At the electroweak scale, a noncanonical 
kinetic term for light neutrinos is generated, which, after canonical normalization, results in 
a nonunitary relation between the flavor and mass eig enstates N = (l - FF*/2)U, where 
U is a unitary matrix and F = vY y M^ 1 . For \F\ > O(0.1), non-negligible non-unitarity 
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FIG. 1: The RG evolution behavior of the three leptonic mixing angles as functions of the energy 
scale n in the SM. We use Mq = 1 TeV and e = 0.3 as input parameters at the seesaw scale. 
The colored regions correspond to values of the leptonic mixing angles within their la confidence 
intervals in the ND case, while the solid, dashed, and dotted curves denote the same confidence 
intervals in the ID case. 



effects could be visible in the near detector of a future neutrino factory, in particular in the 
v n ~^ v t channel, and there are further constraints coming from the universality tests of weak 
interactions, rare leptonic decays, the invisible Z width, and neutrino oscillation data Yl \. 
Thus, in what follows, we restrict ourselves to \Y V \ < 0.3, ensuring full compatibility with 
the nonunitarity constraints. 



A. Single threshold 

For the purpose of illustration, we consider the representative example with M\ ~ M 2 ~ 
M 3 ~ M = 1 TeV, called the single threshold. In Fig. [TJ the RG evolution behavior of 
the three leptonic mixing angles in the SM are shown, in which the input values of Y v at 
/i = M are labeled in each plot. Here, we do not consider the impact stemming from the 
CP- violating phases, but we will comment on that later. The left plot of Fig. [1] shows that, 
in the case of small Yukawa couplings, e.g. \Y„\ < 0.01, there is no visible RG running 
effects on the leptonic mixing angles both below and above the seesaw threshold. This can 
be observed from our analytical formulas, since in the limit x/i — > 0, Q- V j ~ hold for all the 
three leptonic mixing angles. Note that the meaning of e A in the legends of Figs. [TH is that 
the corresponding element is basically negligible compared with elements that are indicated 
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by e. In the middle plot of Fig. [TJ one of the Yukawa couplings in Y v is turned on, i.e., 
yi ~ 0.3, which leads to a significant increase of 6*12 at the GUT scale. On the other hand, if 
Hi is switched off, while yi or y 3 is turned on, as is shown in the right plot of Fig. [TJ # 12 will 
decrease with increasing energy scale due to the opposite sign in front of y\ and 2/2 or y 3 in 
Eq. f)23p . From Fig. [TJ we also find that the RG running of 812 is qualitatively insensitive to 
the hierarchy of the masses of the light neutrinos. This is in agreement with our analytical 
results, since #12 mainly receives corrections from ^ 2 ~ L 1 , which is positive for both possible 
hierarchies. Furthermore, there exist quasifixed points at tt/2 or in the RG running of 812. 
Such a feature can be observed from Eq. fl2"3"j) . where 8\ 2 is proportional to both s 12 and Ci 2 , 
which approaches zero when # 12 ~ or 6 12 ~ vr/2. Finally, due to the lack of a sufficiently 
large enhancement factor, 823 and 8i 3 are relatively stable against radiative corrections. In 
principle, a shift of a few degrees for 623 can be achieved, see, for example, the right plot 
of Fig. [TJ In the limit y2 — 1/3, their contributions to the RG running of 823 and #13 are 
canceled, which can also be seen from Eqs. ( 1241) and ( l2~5lh 

In Fig. El we continue to illustrate the RG running of the leptonic mixing angles in the 
MSSM for both small and large values of tan /3. Compared with the SM, the evolution of 812 
is analogous when the RG running is dominated by Y v . For instance, in the middle and right 
columns of Fig. El #12 approaches gradually its quasi-fixed points in both the small and large 
tan /3 cases. However, if Y v is suppressed, then the leading contributions to the RG running 
of #12 originate from the charged-lepton Yukawa coupling y T , in particular when tan/3 is 
sizable (see, for example, the left column of Fig. ED- This is similar to the RG evolution 
in the effective theory, where the charged-lepton corrections are enhanced by tan/3. Since 
Y e and Y v enter Eq. (I2B1 with different signs, there is a cancellation between them, namely, 
the Y v contributions to the RG running are canceled somewhat by Y e , and vice versa. In 
addition, the RG running of 6*23 and 6*13 may be modified by y T , depending on the specific 
choice of the Yukawa couplings and the value of tan/3. From the lower-right plot in Fig. El 
one can observe that #23 might acquire sizable corrections if tan /3 is large. In order for 813 
to have some visible RG running effects, arrangements of the CP-violating phases have to 
be incorporated, otherwise 8 13 ~ according to Eq. f[24"|) . Therefore, the tiny RG running 
effects in Fig. El come from subleading order corrections. 

The main difference between the RG running of the leptonic mixing angles in the SM 
and the MSSM is that the y T contributions to the running are negligible in the SM, while 
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FIG. 2: The RG evolution behavior of the three leptonic mixing angles as functions of the energy 
scale [i in the MSSM with tan/3 = 10 (upper plots) and tan/3 = 50 (lower plots), respectively. The 
input values of the masses of the heavy singlets and the Yukawa couplings are the same as those 
in Fig. [TJ 



in the MSSM, they are amplified by tan/3. For example, in the case of tan /3 = 50, y T ~ 
m T j[v cos /3) ~ 0.3, indicating a substantial modification of the RG running. In the case 
of tan/3 = 10, one can roughly estimate that y T ~ m T /(wcos/3) ~ 0.1, which leads to a 
relatively small modification. This also reflects the fact that, in the MSSM, the RG running 
of the leptonic mixing angles in the case of a small tan/3 is weaker than that in the case 
of a large tan /3. In addition, the discrepancies between the RG running in the ND and ID 
cases originate from the sign changes in and ^ as shown in Eqs. ( 124)) and (|25|) . For 
this reason, the RG running directions of #23 and #13 are changed if the mass hierarchy of 
the light neutrinos is reversed. Since the beta function of 612 is primarily dominated by 
£12 1 t ne running direction of 6* 12 is independent of the neutrino mass hierarchy, although the 
subleading order corrections may somewhat affect the RG running behavior. 
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FIG. 3: The RG evolution behavior of the three leptonic mixing angles as functions of the energy 
scale // in the SM with the masses of the heavy neutrinos being 1, 10, and 100 TeV, respectively. 
The input values of the Yukawa couplings are the same as those in Fig. [TJ 



B. Multiple thresholds 

In the most general case with a nondegenerate mass spectrum of the heavy singlets, i.e., 
Mi < M 2 < M 3 , the RG running effects between the thresholds may modify the neutrino 
parameters remarkably in the SM. This case is called the multiple thresholds. In Fig. [31 we 
present the typical RG running behavior including three thresholds at 1, 10, and 100 TeV in 
the SM. Note that the chosen values of the three thresholds serve as an order-of-magnitude 
example only. In crossing the thresholds from Mi to M 3 , Q-^i is dramatically suppressed 
no matter the choice of Y v . To identify the distinct threshold corrections, we return to 
Eq. (fTTI) . and observe that the neutrino mass matrix between the seesaw thresholds consists 
of two parts n and Y^M^ 1 Ms(M^)~ 1 Yj . In the SM, the beta functions of these two parts 
have different coefficients in the terms proportional to the gauge couplings and the Higgs 
self-coupling A fl5j |. Therefore, keeping only the gauge coupling and A corrections in the 
corresponding beta functions, after a short-distance running from Mj to Mj_i, the two parts 
of the neutrino mass matrix are rescaled to an and bY u M^ 1 Ms(M^[)~ 1 Yj ', with a ^ b. 
Explicitly, we obtain the mass matrix of the light neutrinos at ji = Mi_i as 



m, 



b [« + YvM^Ms^MD^Yj ] + (a - 6)« 

Mi 



bm, 



+ Ak 

Mi 



Mi 



Mi 

(33) 
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where 

A * 1 (%| + JL,; + a) In . (34) 

16tt 2 V2 2 HT 1 J \ Mi J K J 

One may treat the second term in Eq. fl33l) as a perturbation, arising from the RG running 
between the different thresholds. In contrast to the single threshold case, in which the 
RG corrections to the leptonic mixing angles are due to the large Yukawa couplings Y e 
and Y u , the threshold corrections in this case are related to the gauge couplings and A. 
Their nontrivial flavor structure is the reason for the observed RG running between the 
thresholds. Furthermore, if the masses of the light neutrinos are nearly degenerate, a short- 
distance running will result in sizable corrections to the corresponding leptonic mixing angle 
due to the enhancement factors C7 ■ 

The dramatic decrease of # 12 at around the first seesaw threshold observed in Fig. [3] 
is worth a further comment. First, the discontinuity of the beta function at the seesaw 
threshold is a direct consequence of the present renormalization scheme, in which the heavy 
singlets are decoupled abruptly. If one would instead use a mass-dependent scheme, the 
decoupling, and also the observed RG running, would be smooth. Although the values of 
the renormalized parameters in the vicinity of the thresholds are scheme dependent, the 
total amount of running is less so. Second, the negative slope of the curve between M\ and 
M 2 is arguably an artifact of the diagonality assumption imposed on the neutrino Yukawa 
coupling matrix. This can clearly be observed from the discussion in Sec. IV Dl and seen in 
the corresponding Fig. [7] where an opposite behavior is observed in the top-down approach 
with a nondiagonal Yukawa coupling matrix. Note that a similar discussion applies also to 
the RG running behavior of #23 observed in Fig. [3j 

Apart from the large variation at the seesaw thresholds, one may expect a further increase 
of 9 12 between M 3 and the GUT scale (cf., the middle plot of Fig. [3]), although it is very 
unlikely for the two large leptonic mixing angles #23 and #12 to unify at the GUT scale, 5 unless 
one increases the values of the Yukawa couplings. On the other hand, a smaller #12 could 
be favored at the GUT scale, in particular a Cabibbo-like angle, i.e., 9c ~ 13°. In addition, 
#23 acquires threshold corrections between the second and third thresholds. Although they 
are milder compared with that of #12, the total amount of running between the thresholds 
can be large if M 3 /M 2 is much larger than what it is in our example in Fig. [3J 

5 Indicating a certain flavor symmetry in the lepton sector at the given energy scale. 
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FIG. 4: The RG evolution behavior of the three leptonic mixing angles as functions of the energy 
scale fi in the MSSM with the masses of the heavy neutrinos being 1, 10, and 100 TeV, respectively. 
Similar to Fig. [3l we use tan/3 = 10 for the upper plots and tan/3 = 50 for the lower plots. The 
input values of the Yukawa couplings are the same as those in Fig. [TJ 



In the MSSM, as we have shown in Sec. IHIt there is no mismatch between the RG running 
in the full and the effective theories. Therefore, the running behavior of all three leptonic 
mixing angles should be analogous to the ones in the single threshold case. The numerical 
results are presented in Fig. HJ which are in agreement with our expectations. 

In conclusion, the RG evolutions of the leptonic mixing angles between thresholds may 
lead to sizable corrections to the RG running in the SM, owing to vertex renormalization. 
In general, one should not simply use the matching condition at a common scale, unless the 
mass spectrum of the heavy singlets is considerably degenerate or the mass spectrum of the 
light neutrinos is very hierarchical (i.e., the NH or IH cases). Furthermore, we would like 
to comment on a phenomenologically meaningful situation, in which two of the masses of 
the heavy singlets are close to each other, while the other one is not, e.g., Mi ~ Mi < M3 
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FIG. 5: The RG evolution behavior of the masses of the light neutrinos (left plot) and the CP- 
violating phases (right plot) as functions of the energy scale \x in the ND case. We use Y u = 
diag(0.3, 0.1, 0.1) and Mq = 1 TeV as an example. In the left plot, the solid, dashed, and dotted 
curves correspond to the cases of the SM, the MSSM with tan ft = 10, and tan ft = 50, respectively. 
In the right plot, we show the case of the SM only. 



or Mi < M 2 — M 3 . 6 In such scenarios, the threshold corrections to 812 and #13 are similar 
to those in the general three threshold scenario. However, in the latter case, 823 is free 
of threshold effects, since it only receives visible threshold corrections in the RG running 
between M2 and M3, which can also be seen from Fig. [3j 



C. Neutrino masses and CP-violating phases 

For completeness, we perform the RG running of the neutrino masses and the CP-violating 
phases, despite the fact that there is still a lack of information on leptonic CP violation. As 
shown in Fig. [5] for the ND case, the RG running provides a common rescaling of the masses 
of the light neutrinos, and there is no sudden change along the running direction. It can 
also be seen from Eqs. (!29|) - (|3T1) that the RG running of the masses of the light neutrinos is 
mainly governed by the flavor blind part a v , while the Y v contributions are relatively small. 
There exits a maximal value of the masses of the light neutrinos in the MSSM at the energy 
scale \i ~ 10 10 GeV due to the cancellation between gauge couplings and Yukawa couplings 

6 Such arrangements of the masses of the heavy singlets may come from certain underlying flavor symmetries 
under which the two heavy singlets transform as a doublet. 
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FIG. 6: An example of the RG evolution of the leptonic mixing angles with a nondiagonal Y v in 
the SM. We focus on the NH case and use m\ = 0.01 eV as a sample value. 



in a UJ while such a local maximum does not appear in the SM. Note that similar RG running 
features of the masses of the light neutrinos also exist in the ID case. 

In the single threshold case, the RG running of the CP- violating phases is also dominated 
by the corresponding Yukawa couplings. According to Eqs. (I26 i) -( l28l) . there are enhancement 
factors for both the Dirac and Majorana phases. Interestingly, even if one has a vanishing 
Dirac phase 5 at a certain high-energy scale, a nonvanishing 5 can be generated radiatively 
via the Majorana phases. This can be observed from the right plot of Fig. [5l where 5 = is 
obtained at the GUT scale. As we have pointed out, there is a subtlety at #13 = 0, since the 
definition of 5 loses its meaning. Hence, we choose S13 = 0.01 in our realistic calculations. 
In addition, the Majorana phases p and a run in the same direction, since the coefficients 
of their beta functions possess identical signs at leading order. A complete analysis of the 
RG running of the CP-violating phases could be interesting and useful for model building. 
However, such a study lies beyond the scope of this work. 



D. Nondiagonal Yukawa couplings 

Until now, our numerical analysis has been based on the assumption that both Y e and 
Y v are diagonal at the GUT scale. For the general case with arbitrary Yukawa couplings, it 
is very difficult to obtain analytical RGEs for the neutrino parameters due to the additional 
nonzero elements in Y v . Thus, in this case, a top-down approach is highly favorable, despite 
the technical challenge of accommodating the low-energy data. 
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Not going into details, we focus on adding further support to the comments made in 
Sec. IVBI on the qualitative features of the behavior of 9u and 6*23 among the seesaw thresh- 
olds in the SM as depicted in Fig. [3J In Fig. El we show a numerical example corresponding 
to the NH case with an off-diagonal neutrino Yukawa coupling matrix of the form 



Y, 



e 



\ 



\ 



(35) 



7 



e 4 

e e A 

Unlike in the diagonal case of Fig. |3l one can find directly that #12 increases across the 
seesaw thresholds and it can change significantly even between M2 and M3. Similarly, in 
Fig. El #23 runs in the opposite direction compared with that in Fig. El and most of the 
running is experienced between M\ and M 2 . Note also that the magnitude of the effects 
in Fig. [6] is smaller than in Fig. |3l since a smaller value of mi is adopted for the sake of a 
better convergence of the top-down approach. Therefore, in Fig. |3l the sharp decline of 6*12 
is a consequence of the assumptions made on Y u rather than a generic feature. 



E. Neutrino mixing patterns and flavor symmetries 

While understanding the flavor puzzle is a fundamental task in particle physics, there 
are various flavor symmetric models constructed at high-energy scales, where the flavor 
symmetry is restored. In principle, any predictions from certain flavor symmetries encounter 
radiative corrections, and therefore, it is essential to take into account the RG running effects 
in determining physical parameters at an observable energy scale. In the effective theory, 
it is well known that the RG running effects are negligibly small in the SM. In the MSSM 
with nearly degenerate mass spectrum of the light neutrinos, 9 12 decreases with increasing 
energy scale. In view of these features in the effective theory, there is no way to either unify 
#12 and #23 as suggested in the bimaximal mixing pattern (s 12 = s 2 3 = l/\/2 and Si 3 = 
in the standard parametrization), or arrange them to the tri-bimaximal mixing pattern 
(S12 = l/vo, S23 = 1/V2, and S13 = in the standard parametrization) at high-energy 
scales. Nevertheless, once the RG running effects above the seesaw scale are included, it 
makes sense to realize certain interesting mixing patterns at a unification scale. In Fig. [71 
we illustrate the possibility of realizing the bimaximal and tri-bimaximal leptonic mixing 
patterns at the GUT scale via fine-tuning of the Yukawa couplings in the RGEs. For 
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FIG. 7: Realizations of the bimaximal (left plot) and the tri-bimaximal (right plot) neutrino mixing 
patterns at the GUT scale. Here we only consider the ND case. In addition, for simplicity, we 
use Mo = 1 TeV for all the heavy singlets and take all the CP-violating phases to be zero. The 
solid, dashed, and dotted curves represent the SM, the MSSM with tan/3 = 10 and tan/3 = 50, 
respectively. 



simplicity, we only show the single seesaw threshold case. As examples, in the left plot 
of Fig. [7J we use 

{yui,yu2,yu3) = (0.083, 0.010,0.010), (0.085, 0.040, 0.010), (0.266,0.260,0.010) 

in the SM, in the MSSM with tan /3 = 10, and in the MSSM with tan/3 = 50, respectively, 
whereas in the right plot of Fig. [JJ we use 

(Vvi,yv2,yv3) = (0.034,0.010,0.010), (0.066,0.040,0.010), (0.261,0.260,0.010) 

in the SM, in the MSSM with tan/3 = 10, and in the MSSM with tan/3 = 50, respectively. 
In all cases, flavor symmetric mixing patterns are obtained according to the RG running. 
Thus various mixing patterns can be naturally achieved by adjusting the Yukawa couplings 
in the inverse seesaw model. 



VI. COMPARISON WITH OTHER SEESAW MODELS 

As discussed in Sec. HI there is no difference between the RG evolutions in low-scale 
effective theories. If we stick to the conventional type-I, -II, -III, and inverse seesaw models, 
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at energy scales above the seesaw threshold, the RG running of the neutrino mass matrix 
can be uniformly described by 

m v = N u m u + m v Nl + a u m u , (36) 

where 

N v = C e Y e Yj + C A Y A Yi, (37) 

with Ya being the Yukawa couplings between heavy seesaw particles and lepton doublets. 
The coefficients a v , C e , and C A depend on the specific model. Apparently, sizable Yukawa 
couplings in Eq. fl37|) could give birth to visible RG running effects. In addition, threshold 
effects may induce significant corrections. 

Let us now investigate these two points in detail. In the conventional type-I seesaw model, 
larger Yukawa couplings mean higher energy scale of the right-handed neutrinos, and thus, 
prominent RG running effects could emerge at some super-high-energy scale far above the 
scope of current experiments. If we lower the scale of the right-handed neutrinos by reducing 
the corresponding Yukawa couplings, threshold corrections may play an important role in the 
RG running. However, such a theory is still lacking testability, since the interactions between 
right-handed neutrinos and SM particles are suppressed by the Yukawa couplings, unless we 
use severe fine-tuning or specific assumptions on the model parameters 6j. In the type-II 
seesaw model, one may have sizable Yukawa couplings without facing the problem of lacking 
observability. However, as shown in Ref. \&\, there is no enhancement factor in the type-II 
seesaw framework, and one can hardly have visible RG running effects. For the simplest 
type-II seesaw model, there is only one triplet Higgs, which on the other hand prohibits the 
possibility of threshold corrections. The situation in the type-Ill seesaw model is similar to 
that in the type-I seesaw model, apart from the fact that, due to the charged components, 
there could be visible collider signatures no matter the magnitude of the Yukawa couplings. 

The prominent feature of the RGEs in the inverse seesaw model is that significant RG 
running of the neutrino parameters could occur at lower energy scales without losing the 
testability at the current experiments. The RG running does not spoil the stability of the 
masses of the light neutrinos, but may introduce distinctive corrections to the lepton flavor 
structure. Indeed, this is related to the characteristic of the inverse seesaw model, namely, 
lepton number violation is well separated from lepton flavor violation. 
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Note that, in all cases with visible RG running effects, a nearly degenerate mass spectrum 
of the light neutrinos is required. Otherwise, there is no efficient enhancement factor boosting 
the running. Since the beta functions in Eqs. (l2l?|) - (jHT|) are proportional to the masses of 
the light neutrinos explicitly, a nonvanishing mass cannot be generated via the RG running 
if it is zero at a certain energy scale. 

VII. SUMMARY AND CONCLUSIONS 

In this work, we have performed both analytical and numerical analyses of the RG running 
of the neutrino parameters in the inverse seesaw model. We have shown that, due to the 
sizable Yukawa couplings between light neutrinos and heavy gauge singlets in the inverse 
seesaw model, substantial RG running effects can be naturally obtained even at low-energy 
scales. Such a running distinguishes the inverse seesaw model from other simple seesaw 
models, and may be experimentally tested in near-future experiments. Concretely, we have 
derived very compact analytical RGEs for the neutrino parameters above the seesaw scale. 
Furthermore, a detailed numerical study of the RG running effects on the leptonic mixing 
angles has been carried out in both the SM and the MSSM. In general, we have found 
that there may be significant RG running effects on # 12 and #23, and in particular on # 12 , 
if the mass spectrum of the light neutrinos is nearly degenerate. Furthermore, the running 
between the seesaw thresholds corresponding to a certain hierarchy in the masses of the heavy 
neutrinos can be strong. The RG running effects of a nondiagonal Yukawa coupling matrix 
have also been briefly discussed. We have demonstrated that some phenomenologically and 
theoretically interesting leptonic mixing patterns, the bimaximal and tri-bimaximal patterns, 
can be achieved at a high-energy scale once the RG running effects are taken into account. 
In addition, the RG evolution of neutrino masses and CP-violating phases has been studied 
qualitatively. We have found that the Majorana phases run in the same direction and that 
a nonzero Dirac phase at low-energy scales can be generated from a vanishing one at some 
high-energy scale. 
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Appendix A: Complete RGEs for neutrino parameters 



The complete one-loop RG running of the masses of the light neutrinos are given by 
rhi = 2m 1 {c\ 2 [s 2 13 (s 2 23 yl + c 2 23 yl) + c\ 3 y 2 ] + 2s 12 ci 2 .S23C23.si3C ( 5 (yl - yl) 

+ s2 12 (4 3 Vl + 4 3 vl) + a »} , ( A1 ) 

m 2 = 2m 2 [c\ 2 (c 2 23 y 2 2 + s 2 23 yl) + 2si 2 ci 2 s 23 c 23 s 13 c& (yl - y 2 ) + s 2 12 cl 3 y 2 

+«i 2 s ?3 (4 3 vl + c 2 23 yl) + a v ] , (A2) 
m 3 = 2m 3 { [c\ 3 {s 2 23 yl + c 2 23 yl) + s 2 13 yf] + a v ) , (A3) 

where we have defined % = 167r 2 /i^-, c<5 = cos 5, s$- p = sin(5 — p), and so on. Similarly, 
the full analytical results for the leptonic mixing angles and the Dirac CP-violating phase 
read 

#12 = C2iV-<7 {S23C23S13 (2c s s\ 2 Sp- a - ss + p-a) (yl - yl) 

+ s 12 c 12 s p . a [c\ 3 yl + (s 2 23 s 2 13 - c 2 23 ) y 2 2 - (s 2 23 - c 2 23 s\ 3 ) y 2 3 ] } 
+ C31S12S13S5+P {S12S23C23SP (yl - yl) + ci 2 si 3 s 5+p [(yl - yl) s 2 23 + y\- yl] } 
+ (32Si 3 S6+a {s 2 3C 2 3S CT (yl - yl) c\ 2 - s 12 c 12 s 13 s 8+a [(yl - yl) s 23 + y\ - y 2 3 ] } 

+ -^Cp-a {-523C23-513 (cs+p-a ~ 2c & C p ^ a s\ 2 ) - yfj 

S21 

+ Si2Ci 2 c p „ CT [c\ 3 yl + (s 2 23 s\ 3 - c 2 23 ) y 2 2 - (s 2 23 - c 2 23 s\ 3 ) yf] } 
+ —C5+ P s 12 s 13 {S12S23C23CP (yl - yl) + Ci 2 si 3 c 5+P [(yl - yl) s 2 23 + y{- yl] } 

S>31 

+ -—c& +a s 13 {ci 2 s 2 3C23C CT (yl - yl) - s 12 c 12 s 13 c 5+a [(yl - yl) s 2 23 + y\ - y 2 3 ] } ,(A4) 

S32 

#23 = C3isi 2 s p {c 23 si2S23Sp (vl - vl) - c 12 s 13 s 5+P [(yl - yl) s 2 23 + y\ - yf\ } 
+ (wSa {c 2 12 s 23 c 23 s a (yl - y 2 3 ) + s 12 c 12 s 13 s 5+a [(yl - yl) s 2 23 +yf- yl] } 
+ -^c p si2 {si 2 s 23 C23Cp (yl - yl) - c 12 s 13 c 5+p [(yl - yl) s 2 23 + yf- yl] } 

S31 

+ ^c CT {ci 2 S2 3 C23C ff (yl - yf) + s 12 c 12 s 13 c 5+a [(yl - yl) s 2 23 + y\ - y 2 3 ] } , (A5) 

S32 

#13 = C3lCl 3 S5+p {ci2Si 3 S<5 + p [(vl ~ vl) S l 3 + Vl ~ vl] + Si2C12.S23C23.Sp (yl ~ yl) } 

+ C32S12C13 s s+ a {ci2S 23 C23S CT (yl - yl) + s 12 s 13 s s+a [(yl - yl) sj 3 + yf- yl] } 
+ j-c 13 c s+p {c\ 2 s 13 c 5+ p [(yl - yl) s 2 23 + y\- yl] + si 2 ci 2 S23C 2 3Cp (yl - yl) } 

(,31 

+ -^si 2 ci3C5+ CT {ci 2 S2 3 C23C CT (yl - yl) + s 12 s 13 c 5+a [(yl - yl) s 2 23 + yf- yl] } ,(A6) 

<532 
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{2c p _ ff Si 3 S23 (ss+p~a ~ 2c 5 s 2 12 s p _ a ) (yl - yf) c 2 23 



2Ci2C 23 Sl2 

i T22 /22 2 \ 2 , / 2 22X 21 1 

+ C 12 Si2S2p-2a L _ C 132/l ~~ l S 13 S 23 ~~ C 23j ^2 + \ S 23 ~ C 12 S 23j Vz J C 23 / 

+ 7j — {C12C23S13S23 [((«23 ( s 2«5+2pS? 3 + 2s 2p + S 2 8+2p) ~ S 2p ) yl 

•^Ci2C23Si3S23 

+ (s25+2pSi 3 + 2C 5 S 5+ 2 P - S23 (s 2 5+2pSi 3 + 2s 2 p + S 2 5+2p)) yf) s\ 2 

- (*12 («13 + 1) - 1) S 2<5+2p2/? - S 25 +2p ((vl ~ vl) ^23 + vl)] 

+ s i2 [4a ((^c P ss +P - S12 (s« + 3s 5+2 p)) {yl - yl) s 23 

+ {{2s\ 2 s s+ 2p - 2c p sa +(9 ) yl + 2 (3c p s s+p - s\ 2 (s s + 2s 5+2 p)) yf) 

s 23 ~ 2c 12 C p 

(24. - 1) 

ss+pUi — 2c\ 2 c p ss+ p y1) ~ 2c^2C2 3 C5 + pS2 3 Sp (1/2 — 2/3)] } 

+ ~ — {C12C23S12S13S23 [((2s 2 a + S 2 S+2a) (yl ~ yl) -§23 + S 2a (yl ~ yfj 

•^Ci2C2 3 Si2Sl 3 S23 

+ ( S 13 + l ) S 2S+2aVl + S 2 13 S 2 5+2a (yl ~ yl) S 23 ~ S i3 S 25+2aVl - $28+2^1) 
42 + (2^23 - 1) S2a (yl ~ JJ3) ~ S 2 13 S 2S +2a ((yl ~ yl) ^23 + v\ ~ yf)] 

~ c i2 [( S is (( s s + 3s 5+2ct ) (yl - yl) s\ 3 + (2 (s s + 2s s+2(J ) yl - 2s s+2(J yl) s 2 23 

- 2c a (2sj 3 - l) ss+aljl - 2c ff S5+a2/ 3 ) 

— 2c2 3 C5 +(T S23 s o- (yl — yf)) s? 2 + 2c\ 3 C5 +a s\ 3 s\ 3 s u (y 2 — yl)] } 
+ 7; T~ {~ 2 ( c s+ P -a — 2csc p ^ a s 2 12 ) s 13 s 23 s p ^ a (yl — yl) c\ 3 

^Ci2C2 3 Si2C,21 

- Ci 2 Si2S 2 p_2 CT [~c\ 3 yl - (s 2 13 s 2 23 - c 2 23 ) yl + (4, - c\ 2 s 2 23 ) yl] C23] 

+ ~ j- {C12C23S13S23 [(- (2s 2 p + s 2 «5+2p) (yl - yl) s 2 23 

^Ci2C2 3 Si 3 S2 3 C, 3 i 

+s 2 13 s 2S+ 2p (yl - yl) sj 3 + s 2p yl - 2c s s s+ 2 P yl - s 2 l3 s28 + 2pyl) s i2 

-.2 , „.2 
'23 



+ ( s 12 ( s 13 + !) - l) s 25+2 P yl + S 2 8+2p ((yl ~ yl) s\ 

+S12 [(- ((ss ~ 3s 5+2 p) s\ 2 + Ac s+P s p ) (yl - y 2 3 ) s\ 3 

+ (2 (c 5+p s p - s\ 2 s 5+2p ) yl-2 ((s s - 2s 5+ 2 P ) s\ 2 + 3c,5 + pSp) yl) sl 3 

+ 2c^ 2 C5 + p (2s2 3 — l) s p y\ + 2cl 2 cs +p s p yl) sl 3 + 2c^ 2 C2 3 CpS2 3 s,5 + p (yl — yl)] } 



1 



s 



23 



-2 (5,5+2^2 + (a« - 2s,5 + 2a) yl) s 2 23 - 2c 5+CT (2s2 3 - l) s a yf - 2c s+(T s a yl) 
2c 23 c a s 23 ss+ a (2/2 — %)) s i2 + 2c2 3 c cr s 13 s 23 s ( 5 +(J (y 2 — J/3))] c 12 

+C23S12S13S23 [(S? 3 S 25+ 2 CT (j/a - Vl) S 2 23 + (2s 2ff + S 2S +2a) (vl ~ yl) ^33 
- (s\ 3 + l) S 2 S+2ayl ~ S 2 aVl + s 2 CT y3 + •513^25+2^^3 + S2<5+2 CT 2/f) 
S L - (2^23 - 1) S2o (yl - yl) + s\ 3 S25+2a ( (yl ~ yf) s\ 3 + yf - yl)] C12} . 
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In addition, the analytical RGEs for p and a can be obtained by combining Eqs. (Ilip . (fi~9l) . 
and (|20|) . However, the corresponding formulas are rather lengthy, and therefore, we do not 
list these tedious results here, since the evolution behavior is well described by Eqs. (12"oT) - (|2"S1) . 
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